Relativistic configuration-interaction calculation of energy levels of core-excited states in 

lithium-like ions: argon through krypton 
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Large-scale relativistic configuration-interaction calculation of energy levels of core-excited states of lithium- 
like ions is presented. Quantum electrodynamic, nuclear recoil, and frequency-dependent Breit corrections are 
included in the calculation. The approach is consistently applied for calculating all n = 2 core-excited states 
for all lithium-like ions starting from argon (Z = 18) and ending with krypton (Z — 36). The results obtained 
are supplemented with systematical estimations of calculation errors and omitted effects. 
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I. INTRODUCTION 

Spectroscopy of lithium-like ions has received remarkable 
attention over the last two decades. To a large extent, this at- 
tention was triggered by the famous experiment by Schweppe 
et al. UJ, which demonstrated that the ls^2p-ls^2s energy 
difference in lithium-like ions can be measured with an accu- 
racy of better than 0.1% of the QED contribution. This result, 
recapitulated and surpassed in the following experiments (|2- 
1^, to name a few), made lithium-like ions arguably the best 
playground for testing QED effects in the region of strong nu- 
clear binding field. 

A large number of sophisticated ab initio QED calcula- 
tions have been accomplished during the past years in order 
to advance theory up to the experimental level of accuracy. In 
particular, the two-electron self-energy, two-electron vacuum- 
polarization, and the two-photon exchange corrections were 
calculated |0-[iltl- As a result, the ls^2pj-ls^2s transition 
energies in lithium-like ions are presently among the most 
precise in the medium- and high-Z region. This is not the 
case, however, for the higher excited states of lithium-like 
ions, which received much less attention from theorists. In 
the present work, we aim to correct this drawback. 

The object of our present study are the energy levels of 
the n — 2 core-excited states of medium-Z lithium-like ions. 
Their precise knowledge is required for interpretation of astro- 
physical spectra, as these states contribute to the most promi- 
nent K-shell emission lines, observable in spectra of nearly all 
classes of cosmic X-ray sources (see, e.g., review ifTill ). An- 
other important motivation for studying the K-shell emission 
lfT5i[l6.1 is that it is used for the diagnostics of hot laboratory 
plasmas, particularly those produced in magnetically-confined 
fusion research. 

Accurate theoretical description of transitions involving the 
core-excited states is complicated by two main issues. The 
first one is a large contribution of the QED effects. Indeed, the 
QED effects are strongest for the K-shell electrons and, there- 
fore, for the K-shell transitions. The second problem is that 
the core-excited states are mostly autoionizing states. This 
means that one might expect a strong mixing of the reference 
state with the closely-lying continuum of single excited states 
(i.e., a closed core + an electron in continuum). This inter- 



action with continuum is very difficult to treat accurately in 
theoretical calculations. 

First accurate theoretical data on the core-excited states of 
lithium-like ions were obtained by Vainstein and Safronova 
within the 1/Z-expansion method iTt*]. Later, these states 
were addressed within the multiconfigurational Dirac-Fock 
approach lfl8i - l20tl and, more recently, within the relativistic 
many -body perturbation theory (iT]. In the present work, 
we employ the configuration-interaction (CI) method to ob- 
tain accurate values for the Dirac-Coulomb-Breit energies, 
including the nuclear recoil contribution and the frequency- 
dependent Breit correction. The CI energies are supplemented 
with the QED correction calculated separately in the one- 
electron approximation with a local screening potential. This 
approach is much simpler (and, consequently, less accurate) 
than the full-scale QED calculations reported in Refs. i9l- [l3ll 
for the ls^2pj-ls^2s transitions, but it reproduces the results 
of those more complete calculations remarkably well. Our 
present approach is similar to the one applied previously by 
Chen et al. |22] for calculating the lowest lying n — 2 energy 
levels of lithium-like ions. 

The developed approach is applied for a systematic calcu- 
lation of energy levels of all n = 2 core-excited states for all 
lithium-like ions starting from argon (Z = 18) and ending 
with krypton {Z — 36). For some of these ions, no data on the 
core-excited states were previously available in the literature. 

An important feature of the present investigation is a sys- 
tematic estimation of uncertainties of the calculated results. 
For each particular state and each nuclear charge, we perform 
our CI calculations with a large number (typically, 20-30) of 
different sets of configuration-state functions. By analyzing 
the successive increments of the results obtained with the set 
of configuration-state functions, which is increased in all pos- 
sible directions, we obtain a realistic estimate of how well 
our CI results have converged. Beside the uncertainty of the 
Dirac-Coulomb-Breit energy, we also estimate the uncertainty 
due to the omitted higher-order QED effects. This is done by 
performing the QED calculations with three different screen- 
ing potentials and analyzing the dependence of the results on 
the choice of the potential. 

The remaining paper is organized as follows. In the fol- 
lowing section, we give a brief outline of our implementation 
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of the CI method. The evaluation of the QED correction is 
discussed in Sec.|III] Sec. |IV] presents the numerical details 
of the calculations. Finally, in Sec. [V] we report the results 
of our calculations and compare them with data available in 
the literature. Relativistic units h = c — I and charge units 
e^/AiT = a are used throughout this paper. 



II. CONFIGURATION INTERACTION METHOD 
A. DCB Hamiltonian 

Relativistic Dirac-Coulomb-Breit (DCB) Hamiltonian of an 
iV-electron atom with an infinitely heavy nucleus is given by 



(1) 



where indices i,j — 1, . . . ,N numerate the electrons, hjj is 
the one-particle Dirac Hamiltonian, 



hoii) = a, ■ Pi + {13 ~ l)m + Vnudn) , 



(2) 



a and (3 are the Dirac matrices, p is the momentum operator, 
VJiuc is the binding potential of the nucleus, Vc and Vb are 
the Coulomb and the (frequency-independent) Breit parts of 
the electron-electron interaction. 



Vci^,J) 



2r, 



(3) 



[a^ ■ aj + {a^ ■ r-ij) {otj ■ nj)] , (4) 



with Tij = \ri — rj \ and f = r/r. The operator Hucb acts in 
the space of the positive-energy eigenfunctions of some one- 
particle Dirac Hamiltonian hn, which does not necessarily 
coincides with (|2]i. It is usually convenient to include a part 
of the electron-electron interaction effect already in the one- 
electron orbitals, by introducing a screening potential 14ci (,t) 
and defining ho as 

haii) = Oil ■ pi + {13 - l)m + Kuc(ri) + ^scr(ri) . (5) 

The iV-electron wave function of the system with parity P, 
angular momentum quantum number J, and its projection M 
is represented as a linear combination of configuration-state 
functions (CSFs), 



A^CSF 



^{PJM) = > Cr^{-frPJM) 



r=l 



(6) 



where 7^ denotes the set of additional quantum numbers that 
determine the CSF. The CSFs are constructed as antisym- 
metrized products of one-electron orbitals of the form 



1 f G„(r)xK„m„(r) 
r I zF„(r)x-K„m„(r) 



(7) 



where Xnm the spin-angular spinor 112411 . k — 
(— + 1/2) is the relativistic angular param- 
eter, and m is the angular momentum projection. In the 



present work, we chose the one-electron orbitals i^n to be 
the (positive-energy) eigenfunctions of the one-electron 
Dirac Hamiltonian (|5]l with the screening potential being 
the frozen-core Dirac-Fock potential Vdf, defined by by its 
action on a wave function as 



VDF{ri)i^{ri) 



V / dr2V^+(r2) — 



(8) 



where the index c runs over the core orbitals. 

In the CI method, the energy levels of the system and the 
mixing coefficients in Eq. (|6]l are obtained by solving the 
secular equation 

det{(7,PJAf |i?DCBl7s-PJM) - ,5,4 = . (9) 

The matrix elements of the Hamiltonian between the CSFs 
can be represented as linear combinations of the one- and two- 
particle radial integrals ll25ll . 

{-1rPJM\ Ht,cb\1sPJM) = J2 '^rsiab) I{ah) 

ah 

+ a^^v'^Jl>{abcd)Rk{abcd) . (10) 

k abed 

Here, a, b, c, and d specify the one-electron orbitals, drs and 
Wrs' are the angular coefficients, I{ab) are the one-electron ra- 
dial integrals, and Rk {abed) are the relativistic generalization 
of Slater radial integrals. Further details of our implementa- 
tion of the CI method can be found in our previous study ll26ll . 



B. Corrections to DCB energies 

In the present work, we include two corrections to the DCB 
Hamiltonian ([T]i, namely the nuclear recoil and the frequency- 
dependent Breit corrections. The nuclear recoil correction is 
small and will be treated nonrelativistically. It is conveniently 
separated into two parts, the normal mass shift (NMS) and the 
specific mass shift (SMS). The NMS part can be easily factor- 
ized out and accounted for by multiplying the eigenvalues of 
the DCB Hamiltonian ([T]| by the reduced mass prefactor 



-E^DCB — -Edge , 
m 



(11) 



where /i/m = 1/(1 + m/M) and m and M are masses of the 
electron and the nucleus, respectively. The SMS part of the 
recoil effect is accounted for by adding an additional term to 
the DCB Hamiltonian, 



5Hi 



SMS 



(12) 



The frequency-dependent Breit correction comes from 
the QED theory. It can be obtained by substituting the 
Coulomb Vc and the standard Breit Vb interactions by the 
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full QED electron-electron interaction operator, when calcu- 
lating the matrix elements of the DCB Hamiltonian ([TJ with 
the reference-state wave function(s). Note that in our work, 
contrary to the approach sometimes used in the literature (e.g., 
in Ref. |22]), we do not construct the whole DCB Hamilto- 
nian matrix with the frequency-dependent Breit interaction, 
but apply the frequency-dependent Breit interaction for the 
reference CSFs only. The reason is that the DCB Hamiltonian 
([T]i was derived ll23ll for the standard (frequency-independent) 
Breit interaction only, and inclusion of the frequency depen- 
dency for highly excited CSFs can lead to sizeable spurious 
effects. 

The QED electron-electron interaction operator I{uj) de- 
pends on the energy of the virtual photon uj and has different 
forms in different gauges. In the Feynman gauge, it is given 
by 



-?^Foyn(w, Tij) = a (l - • Ctj) - 

whereas in the Coulomb gauge, it reads 

-a\l+ ("'■'^OK-'^J 



(13) 



(14) 



General QED expressions for the one-photon exchange cor- 
rection for a given reference state Ii27i1 involve two kind of 
matrix elements, the direct ones of the form {ab\I\ab) and the 
exchange ones of the form {ba\I\ab). The energy of the virtual 
photon is fixed by w = in the direct terms and by a; — Sa—Sb 
in the exchange terms. In the case of several equivalent ref- 
erence CSFs, there are also non-diagonal terms and conse- 
quently, matrix elements of a general form {ab\I\cd). For 
them, we used an energy symmetrized expression i28ll 



1 



{ab\I\cd) -> - [{ab\Iiea 



.)\cd) 



{ab\I{eb - £d)\cd)] . 

(15) 

It can be easily verified that, for a single reference state, the 
matrix elements of the operator / are gauge independent. In 
the case of several reference-state configurations, gauge in- 
variance is not exact, but we checked that the residual is com- 
pletely negligible in each particular case. 



III. RADIATIVE QED CORRECTIONS 

In this section, we deal with corrections to the DCB ener- 
gies due to QED effects. For the purpose of the present in- 
vestigation, it will be sufficient to take into account only the 
dominant part of the radiative QED effects, given by the one- 
electron self-energy and the vacuum polarization calculated in 
a local screening potential. 

The total QED contribution for a given many-electron ref- 
erence state is obtained by summing the QED contributions 



from the one-electron orbitals, weighted by their fractional oc- 
cupation numbers as obtained from the eigenvectors of the CI 
calculation, 

SEqed ^J2'i<^ [(alI]sE(£a)la) + {a\Vvp\a)] , (16) 



where index a runs over all one-electron orbitals contributing 
to the given many-electron state, qa is the occupation number 
of the one-electron orbital, Egg is the self-energy operator, 
is the Dirac energy of the one-electron state a, and Vyp is the 
vacuum polarization potential. 

The one-electron self-energy correction in a general lo- 
cal potential is calculated rigorously, to all orders in the nu- 
clear binding strength parameter Za. The method used in 
the present work is quite close to the one developed by Blun- 
dell and Snyderman f?^]. We, however, introduced several 
improvements to their original calculation scheme. First, we 
used the dually kinetically balanced i?-spline basis set i30ll . 
which yielded significantly better con verg ence than the origi- 
nal S-spline method by Johnson et al. II31I1 . Second, we evalu- 
ated the many-potential electron propagator G'^^-' by forming 
the difference 



G(2+) = G- G(") - G(i) , 



(17) 



where G, G^"), and G^^) are the full electron propagator, 
the free propagator, and the one-potential propagator, respec- 
tively. The one-potential propagator was calculated by taking 
a numerical derivative of the B-spline basis set with respect 
to the nuclear charge Z, 



(18) 



We found that this approach yields results equivalent to those 
obtained by the direct (triple) summation over the spectrum 
but is much faster in actual computations. 

The one-electron vacuum-polarization contribution was 
straightforwardly calculated as an expectation value of the 
Uehling and Wichmann-Kroll potentials with the one-electron 
orbitals. The Wichmann-Kroll potential was evaluated by the 
approximate formulas tabulated in Ref. | l32ll . 

In our calculations of the one-electron self-energy and 
vacuum-polarization corrections, we assume that the initial 
state and all intermediate states in the electron propagator are 
the eigenstates of the Dirac Hamiltonian Q with a screening 
potential Vscr- The screening potential allows us to account 
for the dominant part of the electron-electron interaction and 
should be chosen with care. It would be natural to use the 
same screening potential as in the CI part of our calculations, 
but this is not possible since the Dirac -Fock potential ^ is 
non-local. Because of this, we introduce three different local 
screening potentials, whose eigenstates and eigenvalues are 
quite close to the Dirack-Fock ones. 

Each potential is constructed into two steps. First, we solve 
the Dirac-Fock problem (for the center-of-gravity of all equiv- 
alent relativistic configurations) and find the Dirac-Fock wave 
functions. Second, we use these wave function to set up our 
local screening potentials. 
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The simplest choice is the potential induced by the charge 
density of all other electron orbitals except the reference one. 
This is the core-Hartree (CH) potential defined, for a given 
one-electron state a, as 



dr'- 



1 



7 -^Pa{r') , 

(r, r') 



(19) 



where pa is the density of all one-electron orbitals excluding 



'Fbir)] , (20) 



with Gb and Fi, being the upper and the lower components of 
the Dirac-Fock solution. Note that the CH potential defined 
in this way is different from the one frequently encountered in 
the literature because we do not require self-consistency for 

VcH,a- 

The second choice of the screening potential is a variant of 
the Kohn-Sham potential SH, defined as 



VKsir) 



a I dr' )- 

/o max(r, r') 



2 a 

3 r 



81 
32^ 



r p{r) 



1/3 



a 
r 



-(Ar) 



where p{r) is the total Dirac-Fock charge density 

p{r)=J2qt [Glir)+F^{r)] . 



(21) 



(22) 



The second term in the right-hand side Eq. ( ITTT i is the ex- 
change correction derived from density functional theory, 
whereas the third term is a kind of the Latter correction Il35[l . 
which restores the physical asymptotical behaviour of the KS 
potential at large distances. The parameter A in the Latter cor- 
rection was chosen to be about Za/10. Again, the KS poten- 
tial is different from the one encountered in the literature 
because we do not require self-consistency in its definition. 

The third screening potential employed in this work is re- 
ferred to as the Localized Dirac-Fock (LDF) potential. For a 
given state a, the LDF potential l^DF,a is obtained by invert- 
ing the Dirac Fock equations with the known solutions Ga and 
Fa and then smoothing out the result in the vicinity of zeros 
of Ga- The LDF potential was first introduced for calculations 
of the QED corrections to the parity-nonconserving transition 
amplitudes in caesium [36., 371. In the present work, we em- 
ploy the variant of the LDF potential described in detail in 
Ref. 

According to our experience and comparison with the re- 
sults of more complete calculations, the LDF potential yields 
better results in calculations of the QED effects than the KS 
and CH potentials. Therefore, we will use this potential 
for obtaining the final results for the QED correction in our 
present calculation. 



TABLE I: Contributions to the Coulomb energy (in a.u.) for the 
ls2s2p'^Py2 state of lithium-like iron {Z = 26), for the infinitely 
heavy nucleus. The values listed after the first row are the increments 
obtained on successively adding configurations while increasing the 
maximal value of the orbital quantum number Lmax and enlarging 
the size of the one-electron basis. Label SD stands for the single 
and double excitations, T denotes the triple excitations, Ua is the 
number of B-splines in the one-electron basis set, e is energy the 
cutoff parameter, see text for the details. 

SE 



SD, Tia = 30, e = 4 

1 

2 

3 

4 

5 

6 

7 

> 8 

SD, na = 40, e = 8 

1 

2 

3 

> 4 

SD, ria = 50, e = 16 

1 

2 

T, = 25, e = 4 

1 

2 

3 

> 4 

Final result 



-497.762 342 48 
-0.008 719 04 
-0.000 662 72 
-0.000 124 13 
-0.000 034 81 
-0.000 012 45 
-0.000 005 18 
-0.000 005 55 

-0.000173 23 
-0.000 00134 
-0.000 001 13 
-0.000 00107 

-0.000 006 46 
-0.000 000 43 

-0.000 000 22 
-0.000 000 77 
-0.000 000 09 
-0.000 000 04 

-497.772 091 (7) 



TABLE II: Contributions to the frequency-independent Breit energy 
(in a.u.) for the ls2s2p^Pj°/2 state of lithium-like iron (Z = 26), 
for the infinitely heavy nucleus. Notations are the same as in Table|I] 



5E 



SD, na = S0,e = 4 

1 

2 

3 

4 

5 

> 6 

SD, na=40,e = 8 

1 

2 

3 

SD, Ua = 50, e = 16 
1 

Final result 



0.065 471 70 
-0.000 330 01 
-0.000 050 23 
-0.000 014 47 
-0.000 005 98 
-0.000 00915 

-0.000 00176 
-0.000 000 43 
-0.000 000 25 

0.000 004 49 
0.065 064 (9) 
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TABLE III: QED contributions (in a.u.) for the Is2s2p Pfyj ~ Is 2s S1/2 transition in lithium-like iron {Z — 26), for different screening 
potentials, q denotes the occupation number of the single electron orbitals as obtained from the CI calculations. SE, VP, and QED label the 
self-energy, the vacuum-polarization, and the total QED contribution, respectively. 









4 po 












4 po 2 c. 


State 


q 


SE 


VP 


QED 


q 


SE 


VP 


QED 





LDF potential: 
Is 1.000 
2s 1.000 
2pi/2 0.943 
2p3/2 0.057 
Sum 

KS potential: 
Is 1.000 
2s 1.000 
2pi/2 0.943 
2p3/2 0.057 
Sum 

CH potential: 
Is 1.000 
2s 1.000 
2pi/2 0.943 
2p3/2 0.057 
Sum 

Final result 



0.15688 
0.01974 
-0.00052 
0.00005 



0.15512 
0.01991 
-0.00053 
0.00005 



0.15674 
0.01956 
-0.00051 
0.00005 



-0.01337 
-0.00158 
-0.00001 
0.00000 



-0.01326 
-0.00159 
-0.00001 
0.00000 



-0.01337 
-0.00157 
-0.00001 
0.00000 



0.14350 
0.01816 
-0.00053 
0.00005 



0.14186 
0.01832 
-0.00054 
0.00005 



0.14338 
0.01799 
-0.00052 
0.00005 



2.000 0.30432 
1.000 0.01894 



2.000 0.30287 
1.000 0.01922 



2.000 0.30384 
1.000 0.01894 



-0.02610 
-0.00153 



-0.02600 
-0.00154 



-0.02607 
-0.00153 



0.27822 
0.01741 



0.27687 
0.01768 



0.27778 
0.01741 



-0.13472 

0.00075 
-0.00053 

0.00005 
-0.13445 

-0.13501 

0.00064 
-0.00054 

0.00005 
-0.13487 

-0.13440 

0.00058 
-0.00052 

0.00005 
-0.13429 

-0.13445 (43) 



IV. NUMERICAL DETAILS 

In the present work, we use our implementation of the CI 
method employed previously in Refs. lES. [38ll for the eval- 
uation of the hyperfine splitting in Li and Be+. The one- 
electron basis was constructed with help of the dual-kinetic- 
balance B-spline method Ii30il . The screening potential in the 
one-electron Hamiltonian (|5]l was taken to be the frozen-core 
DF potential (O, with the Is^ core in the case of the ground 
and valence-excited states and the Is core in the case of the 
core-excited states. The standard Fermi model of the nuclear 
charge distribution was used to represent the nuclear potential. 

For a given number of B-splines Ua, all Dirac eigenstates 
were taken with the energy < e„ < mc^(l + Zae) and 
the orbital quantum number I < Lmax, where the value of the 
energy cutoff parameter e was varied between 4 and 16 and 
-^max, between 1 and 7. In order to check the convergence 
of our CI results with respect to the number of partial waves 
included and the size of the basis, we perform a set of calcu- 
lations with different number of basis functions and analyze 
the successive increments as the size of the basis is enlarged. 
In our analysis, we study the dependence of the results ob- 
tained on the parameters Ua, e, and £max, and the type of vir- 
tual excitations included. In all cases relevant for the present 
study, the contribution of triple excitations was found to be 
very small. We thus perform our calculations first with in- 
cluding single and double excitations only and then adding the 
contribution of the triples calculated separately with a smaller 
basis. The analysis of the convergence of the partial-wave 
expansion is supplemented with an estimation of the omitted 



tail, which is obtained by a polynomial least-squares fitting of 
the increments in 1/L 

An example of our CI calculations is presented in Table U 
for the Dirac-Coulomb energy and in Table for the Breit 
correction. The results reported in these tables are obtained by 
performing calculations with 27 different sets of basis func- 
tions. Each set is characterized by parameters Ua, e, and 
Lmax- E.g., for Ua = 30, 6 = 4, and imax = 7, the one- 
electron basis reads 17s 16p 16d 15/ 15gl5h 14z 14A; , where 
16p means 16 lowest-lying orbitals of the pi/2 symmetry and 
16 lowest-lying orbitals of the P3/2 symmetry, etc. Generating 
all possible single- and double-excited CSF's with this set of 
one-electron orbitals leads to a basis of about A'^csf = 30, 000 
functions. As another example, for Ua = 50, e — 16, and 
-^max = 2, the set of one-electron orbitals is 29s29p28d, 
leading to the basis of about A^csf = 25, 000 functions. 

Analysis of data presented in Tables U and allows us to 
give a reliable estimate of the uncertainty of our CI calcu- 
lations. In the case shown in these tables, the dominant er- 
ror comes from the omitted tail of the partial-wave expan- 
sion. For higher core-excited states, however, the dominant 
error often comes from the convergence of the one-electron 
basis of np and nd symmetries, which is explained by strong 
mixing of the reference (autoionizing, core-excited) states 
with the closely-lying continuum of single-excited Is^np and 
Is^nd virtual states. In some cases (particularly, for the 
states), we had to employ large one-electron basis sets, up to 
40s 39p39d {Nqsp — 47, 000 functions), in order to ensure 
the convergence of CI results. 

We now discuss our calculations of the radiative QED cor- 
rections. TableHnipresents a detailed breakdown of individual 
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QED contributions for the Is2s2p '^P°/2 - ls^2s transi- 
tion in lithium-like iron. The calculation is performed for the 
LDF, KS, and CH potentials. All three results for the transi- 
tion energy agree very well with each other Remarkably, the 
agreement is much better for the energy difference than for 
each energy level separately. The same situation occurs for 
the other transitions from the core excited states to the ground 
state. 

As a final result for the QED contribution, we take the value 
obtained with the LDF potential (since this potential yields 
the best agreement with the results of more complete calcula- 
tions). The uncertainty of the QED contribution was estimated 
as the maximal deviation of the KS and CH values from the 
LDF one. We checked that such estimation of uncertainty is 
consistent with the results of more sophisticated QED calcu- 
lations available for the ls^2pj states iB 111 [Hi ■ The set of 
calculations presented in Table Ullltakes about 30h of the pro- 
cessor time on a modern workstation. 



V. RESULTS 

Table II V [ presents our calculation results for the energy lev- 
els of the ls^2s, ls^2p, and Is2l2l' states of Hthium-like 
ions, from argon (Z = 18) to krypton (Z — 36). The total 
energies are listed in the table for the ground state, whereas 
for all other states, relative energies with respect to the ground 
state are given. Note that for the ls^2pj-ls^2s energy differ- 
ences, much more accurate theoretical results are available in 
the literature Illl fTsll . We nethertheless present our own re- 
sults for the ls^2pj states in the table since they (i) provide 
us with an independent check of accuracy of our calculations 
and (ii) allow us to deduce the transition energies between the 
core-excited and the ls^2pj states. 

The third and fourth columns of Table HV] present our re- 
sults for the Dirac-Coulomb and Breit energies (multiplied by 
the reduced mass prefactor). For each value of the nuclear 
charge Z and each state, the CI calculation was performed 
with 20-30 different basis sets. After that, an extrapolation of 
the partial-wave expansion to L,„ax — oo was performed and 
an estimation of the uncertainty was made by analyzing the 
results as illustrated by Tables HI and HIl The procedure for the 
error estimation is fully automatized and done consistently in 
the same way for all Z and all states. Each set of CI calcula- 
tions for one value of Z and one state takes about 20-50h of 
processor time on a modern workstation. 

The fifth column of Table |IV] presents our calculation re- 
sults for the QED correction. For each value of Z and each 
state, the calculation was performed with three screening po- 
tentials. The results listed in the table are obtained with the 
LDF potential. The maximal deviation of the KS and CH re- 
sults from the LDF one was taken as the uncertainty of the 
QED correction. Small additional contributions to the en- 
ergy levels from the specific mass shift and the frequency- 
dependent Breit correction are listed in the sixth and seventh 
columns of Table II VI respectively. 

One may note that the uncertainties of our calculation re- 
sults presented in Table IIV I have some irregularities along the 



Z sequence. For some values of Z the estimated uncertainty 
turns out to be much larger than for the neighbouring ones. 
This is caused by a non-uniform convergence of the results 
with respect to the size of the basis set, which is due to the 
interaction of the reference autoionizing state with the con- 
tinuum of single-excited states. In the CI calculations, the 
continuum is discretized and, if a continuum state in the basis 
happens to be nearly degenerate in energy with the reference 
state, a sizeable non-physical contribution may arise. (It is 
non-physical since the true Green function G(e) is known to 
be regular at energies e > m.) When this occurs, we observe 
irregularity in the convergence of our results and have to in- 
crease the uncertainty accordingly. Note that for a given state, 
we use exactly the same basis for all values of Z listed in Ta- 
ble |IV] and estimate the uncertainty exactly in the same way, 
without any smoothing of energies or uncertainties along the 
Z sequence. 

Our final theoretical results for the wavelengths of the 22 
strongest \s2l2l' — s- ls^2l' transitions in lithium-like ions are 
summarized in Table IV] The transitions are labelled from "a" 
to "v", following the standard notations as by Gabriel fs^. 

We now turn to the comparison of our calculation results 
for the three most important elements, namely argon, iron, and 
krypton, with data previously reported in the literature. Since 
the amount of the experimental and theoretical data available 
for these ions is rather voluminous, we decided to restrict our- 
selves to a comparison with the relativistic many-body pertur- 
bation calculation by Safronova and Safronova ll2Tll and with 
the extensive compilations of the theoretical and experimen- 
tal world data [40-421. The calculation by Safronova and 
Safronova is apparently the best of the previous theoretical 
studies of the core-excited states and the compilations f4C|- 
|42] are the primary sources for the NIST recommended spec- 
tral data for these ions ll43ll . In the case of argon, the values 
reported in the compilation by Saloman (^\ include the re- 
sults by Safronova and Safronova but in addition provide esti- 
mations of uncertainties, so we compare only with Saloman's 
compilation. 

The comparison for argon, iron, and krypton is presented 
in Tables IVII IVIII and IVIIII respectively. In the case of ar- 
gon, the compilation |40] includes estimations of uncertain- 
ties. The agreement is nearly always within the given error 
bars but our results are significantly more accurate. For iron 
and krypton, there are no estimations of uncertainty given in 
previous works, but the probable errors might be guessed from 
the differences between different results. Agreement of our re- 
sults with results by Safronova and Safronova is in most cases 
better than that with the NIST compilation data, especially in 
the case of krypton. 

Comparison of our calculation results for other ions with 
the NIST spectral data ll43ll looks similar to what can be seen 
in Tables |VH IVIII and IVIIII Root-mean-square deviation of 
energies of core-excited states varies from 0.002 Ry for argon 
to 0.05 Ry for krypton and grows smoothly with the increase 
of the nuclear charge. The only large deviation is found for 
the ls2s2p^P°y2 state of copper, for which NIST database 
reports 615.26 Ry whereas our result is 612.319(1) Ry. This 
deviation is most probably due to a misprint in the NIST 



7 



database. No data on energies of core-excited states is avail- 
able in the database for Z = 30 and Z = 32 — 35 and only 
few levels are available for Z = 20 and 21. 

In summary, we performed a large-scale relativistic CI cal- 
culation of the energy levels of all n = 2 core-excited states 
of lithium-like ions from argon to krypton. The CI results 
include the nuclear recoil contribution and the frequency- 
dependent Breit correction. The QED correction was calcu- 
lated separately in the one-electron approximation with a lo- 
cal screening potential. The calculation results are supple- 
mented with a systematic estimation of uncertainties. The un- 
certainties of the CI values were evaluated by analyzing the 
successive increments of the results obtained with the set of 
configuration-state functions increased in all possible direc- 
tions. The uncertainty of the QED corrections were evalu- 



ated by performing the QED calculations with three different 
screening potentials and analyzing the dependence of the re- 
sults on the choice of the potential. 
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TABLE IV: Energy levels of lithium-like ions from argon through krypton, in Ry- 
dbergs, 1 Ry = 109 737.315 685 39 (55) cm-\ Separately listed are the Dirac- 
Coulomb energy, the frequency-independent Breit correction, the QED correc- 
tion, the specific mass shift (SMS) correction, and the frequency-dependent Breit 
(Breit') correction. The total energies are presented for the ground state, whereas 
for all other states, the energies relative to the ground state are given, n/m is the 
reduced mass prefactor and R is the nuclear root-mean-square charge radius. 



Term 




J 


Coulomb 


Breit 


QED 


SMS 


Breit' 


Total 




Ar(Z = 18), 1 -yu/m 


= 0.00001373, R = 


3.427 fm 
















1/2 


-696.0032 


0.1504 


0.1676 


0.0001 


-0.0001 


-695.6852 


(4) 




2po 


1/2 


2.3412 


0.0113 


-0.0098 


-0.0006 


-0.0003 


2.3417 


(6) 






3/2 


2.5867 


-0.0015 


-0.0092 


-0.0006 


-0.0004 


2.5749 


(8) 






1/2 


226.5266 


-0.1362 


-0.0669 


-0.0001 


0.0001 


226.3234 


(17) 


lsCS)2s2p{^P°) 


4po 


1/2 


227.0887 


-0.1094 


-0.0753 


-0.0007 


-0.0002 


226.9032 


(4) 






3/2 


227.1777 


-0.1359 


-0.0751 


-0.0007 


-0.0003 


226.9658 


(4) 






5/2 


227.3442 


-0.1469 


-0.0747 


-0.0007 


-0.0003 


227.1215 


(3) 




2po 


1/2 


228.9723 


-0.1355 


-0.0751 


0.0003 


0.0004 


228.7624 


(9) 






3/2 


229.0991 


-0.1400 


-0.0748 


0.0005 


0.0005 


228.8852 


(4) 


lsCS)2p\^P) 


4p 


1/2 


229.6458 


-0.1074 


-0.0852 


-0.0013 


-0.0005 


229.4514 


(6) 






3/2 


229.7494 


-0.1074 


-0.0850 


-0.0013 


-0.0005 


229.5553 


(6) 






5/2 


229.8938 


-0.1384 


-0.0847 


-0.0013 


-0.0006 


229.6689 


(6) 


1sCS)2s2pCP'') 


2po 


1/2 


229.8117 


-0.1136 


-0.0749 


-0.0005 


-0.0001 


229.6226 


(39) 






3/2 


229.8798 


-0.1322 


-0.0748 


-0.0006 


-0.0002 


229.6720 


(43) 


lsCS)2p^CD) 




3/2 


230.9327 


-0.1099 


-0.0849 


-0.0006 


-0.0001 


230.7372 


(33) 






5/2 


230.9812 


-0.1507 


-0.0847 


-0.0007 


-0.0003 


230.7449 


(35) 


lsCS)2p\^P) 


2p 


1/2 


231.2057 


-0.1479 


-0.0845 


0.0005 


0.0005 


230.9743 


(12) 






3/2 


231.4382 


-0.1464 


-0.0845 


0.0004 


0.0004 


231.2081 


(7) 


lsCS)2p^CS) 




1/2 


232.5676 


-0.1241 


-0.0823 


-0.0006 


-0.0002 


232.3605 


(133) 



K(Z = 19), 1 - /u/m = 0.00001408, R - 
ls^2s 



ls^2f> 



ls2s^ 



2po 



lsCS)2s2p{''P°) ^P° 

lsCS)2s2p{^P°) 
lsCS)2p\^P) ^P 

ls(25)2s2p(ip°) 



1/2 
1/2 
3/2 
1/2 
1/2 
3/2 
5/2 
1/2 
3/2 
1/2 
3/2 
5/2 
1/2 
3/2 



-777.9530 
2.4910 

2.8025 
253.5172 
254.1116 
254.2222 
254.4336 
256.1274 
256.2908 
256.8569 
256.9900 
257.1686 
257.0338 
257.1170 



= 3.435 fm 
0.1782 
0.0137 

-0.0018 
-0.1614 
-0.1296 
-0.1612 
-0.1742 
-0.1593 
-0.1663 
-0.1273 
-0.1272 
-0.1642 
-0.1357 
-0.1559 



0.2021 
-0.0119 

-0.0112 
-0.0808 
-0.0910 
-0.0907 
-0.0902 
-0.0908 
-0.0904 
-0.1031 
-0.1027 
-0.1024 
-0.0906 
-0.0904 



0.0001 
-0.0007 

-0.0006 
-0.0001 
-0.0008 
-0.0008 
-0.0008 
0.0003 
0.0005 
-0.0015 
-0.0015 
-0.0014 
-0.0005 
-0.0007 



-0.0001 
-0.0004 

-0.0005 
0.0001 
-0.0003 
-0.0003 
-0.0004 
0.0005 
0.0005 
-0.0006 
-0.0006 
-0.0007 
-0.0001 
-0.0002 



-777.5727 
2.4918 

2.7882 
253.2750 
253.8901 
253.9691 
254.1681 
255.8781 
256.0351 
256.6246 
256.7579 
256.8999 
256.8070 
256.8698 



5) 
7) 

7) 
17) 

4) 

4) 

4) 

9) 

4) 

7) 

7) 

7) 

36) 

40) 





2n 

U 




ZOo . ZOZo 


— u.iouy 


— u 


1 097 


— u 


UUUO 


n 

— u 


UUUl 


ZOo 










5/2 


258.3223 


-0.1784 


-0 


1024 


-0 


0008 


-0 


0003 


258 


0404 


(22) 


lsfS)2p2(3p) 


2p 


1/2 


258.5281 


-0.1750 


-0 


1029 





0006 





0006 


258 


2514 


(7) 






3/2 


258.8289 


-0.1726 


-0 


1022 





0005 





0005 


258 


5551 


(8) 


15(^5)2^2(15) 




1/2 


260.0040 


-0.1472 


-0 


1003 


-0 


0007 


-0 


0002 


259 


7557 


(40) 


Ca(Z = 20), 1-fi/m 


= 0.00001373, R = 


3.476 fm 




















ls^2s 




1/2 


-864.5349 


0.2094 





2412 





0001 


-0 


0002 


-864 


0843 


(6) 


ls^2p 




1/2 


2.6419 


0.0163 


-0 


0143 


-0 


0007 


-0 


0004 


2 


6426 


(8) 






3/2 


3.0317 


-0.0021 


-0 


0135 


-0 


0007 


-0 


0006 


3 


0148 


(8) 






1/2 


282.0476 


-0.1895 


-0 


0965 


-0 


0001 





0001 


281 


7615 


(17) 






1/2 


282.6744 


-0.1522 


-0 


1088 


-0 


0008 


-0 


0003 


282 


4123 


(4) 






3/2 


282.8100 


-0.1897 


-0 


1085 


-0 


0008 


-0 


0004 


282 


5106 


(5) 






5/2 


283.0754 


-0.2046 


-0 


1079 


-0 


0008 


-0 


0005 


282 


7617 


(4) 


ls{ S)2s2p{ P ) 


2 po 


1/2 


284.8243 


-0.1855 


-0 


1086 





0003 





0005 


284 


5311 


(9) 






3/2 


285.0324 


-0.1958 


-0 


1081 





0006 





0006 


284 


7297 


(5) 


lsCS)2p^CP) 


4p 


1/2 


285.6153 


-0.1494 


-0 


1233 


-0 


0016 


-0 


0007 


285 


3403 


(8) 






3/2 


285.7844 


-0.1494 


-0 


1229 


-0 


0016 


-0 


0008 


285 


5098 


(8) 






5/2 


286.0022 


-0.1931 


-0 


1225 


-0 


0016 


-0 


0008 


285 


6842 


(8) 


1sCS)2s2pCP°) 


2po 


1/2 


285.8054 


-0.1608 


-0 


1083 


-0 


0005 


-0 


0001 


285 


5357 


(33) 






3/2 


285.9054 


-0.1822 


-0 


1081 


-0 


0008 


-0 


0003 


285 


6141 


(38) 


ls(\S)2p'CD) 




3/2 


287.1271 


— n 1 ^47 


— n 


1228 


— n 




— n 


nnni 




84SS 








5/2 


287.2249 


-0.2093 


-0 


1224 


-0 


0009 


-0 


0004 


286 


8919 


(15) 


lsCS)2p^CP) 


2p 


1/2 


287.4004 


-0.2050 


-0 


1231 





0006 





0007 


287 


0736 


(8) 






3/2 


287.7846 


-0.2015 


-0 


1221 





0005 





0005 


287 


4620 


(8) 


lsCS)2p\^S) 




1/2 


289.0017 


-0.1730 


-0 


1200 


-0 


0007 


-0 


0003 


288 


7077 


(25) 


Sc(Z = 21), 1-n/m-- 


= 0.00001221, R = 


3.544 fm 




















ls^2s 




1/2 


-955.7643 


0.2440 





2853 





0001 


-0 


0002 


-955 


2351 


(6) 


ls'^2p 


2po 


1/2 


2.7937 


0.0190 


-0 


0171 


-0 


0007 


-0 


0005 


2 


7944 


(9) 






3/2 


3.2760 


-0.0026 


-0 


0160 


-0 


0007 


-0 


0007 


3 


2560 


(9) 


Is2s2 




1/2 


312.1222 


-0.2208 


-0 


1142 


-0 


0001 





0001 


311 


7872 


(16) 


lsCS)2s2p{^P°) 


4po 


1/2 


312.7818 


-0.1773 


-0 


1289 


-0 


0008 


-0 


0003 


312 


4744 


(5) 






3/2 


312.9458 


-0.2214 


-0 


1286 


-0 


0008 


-0 


0004 


312 


5946 


(5) 






5/2 


313.2757 


-0.2384 


-0 


1278 


-0 


0008 


-0 


0006 


312 


9081 


(5) 


lsCS)2s2p{^P'') 


2po 


1/2 


315.0672 


-0.2140 


-0 


1287 





0003 





0006 


314 


7254 


(10) 






3/2 


315.3299 


-0.2287 


-0 


1280 





0006 





0007 


314 


9745 


(5) 


ls(25)2p2(3p) 


4p 


1/2 


315.9258 


-0.1741 


-0 


1462 


-0 


0016 


-0 


0008 


315 


6032 


(9) 






3/2 


316.1389 


-0.1740 


-0 


1456 


-0 


0016 


-0 


0009 


315 


8168 


(9) 






5/2 


316.4006 


-0.2254 


-0 


1452 


-0 


0015 


-0 


0010 


316 


0276 


(9) 


1sCS)2s2pCP°) 


2po 


1/2 


316.1327 


-0.1893 


-0 


1282 


-0 


0005 


-0 


0001 


315 


8146 


(32) 






3/2 


316.2512 


-0.2112 


-0 


1280 


-0 


0008 


-0 


0004 


315 


9109 


(38) 


lsCS)2p^CD) 




3/2 


317.5616 


— u. ±o±o 


n 
u 




n 




— n 


nnni 


■^1 7 










5/2 


317.6965 


-0.2434 


-0 


1451 


-0 


0008 


-0 


0005 


317 


3067 


(13) 


lsCS)2p^CP) 


2p 


1/2 


317.8285 


-0.2383 


-0 


1459 





0006 





0008 


317 


4456 


(9) 






3/2 


318.3134 


-0.2333 


-0 


1447 





0005 





0006 


317 


9364 


(9) 


Is(25)2p2(i5) 




1/2 


319.5686 


-0.2019 


-0 


1423 


-0 


0007 


-0 


0003 


319 


2234 


(20) 


Ti (Z = 22), 1 - n/m = 


= 0.00001144, R = 


3.591 fm 




















ls^2s 




1/2 


-1051.6554 


0.2823 





3347 





0001 


-0 


0002 


-1051 


0386 


(7) 


ls^2p 


2po 


1/2 


2.9466 


0.0222 


-0 


0202 


-0 


0007 


-0 


0006 


2 


9473 


(10) 






3/2 


3.5371 


-0.0030 


-0 


0189 


-0 


0007 


-0 


0008 


3 


5136 


(10) 


ls2s^ 




1/2 


343.7454 


-0.2553 


-0 


1341 


-0 


0001 





0001 


343 


3559 


(16) 


1sCS)2s2pCP'') 


ipo 


1/2 


344.4377 


-0.2051 


-0 


1514 


-0 


0009 


-0 


0004 


344 


0800 


(6) 






3/2 


344.6340 


-0.2566 


-0 


1511 


-0 


0008 


-0 


0005 


344 


2250 


(5) 






5/2 


345.0404 


-0.2758 


-0 


1501 


-0 


0008 


-0 


0007 


344 


6130 


(5) 


lsCS)2s2p{^P°) 


2po 


1/2 


346.8601 


-0.2450 


-0 


1512 





0003 





0007 


346 


4649 


(11) 






3/2 


347.1888 


-0.2651 


-0 


1504 





0006 





0008 


346 


7748 


(6) 


lsCS)2p^CP) 


4p 


1/2 


347.7928 


-0.2014 


-0 


1716 


-0 


0016 


-0 


0009 


347 


4173 


(10) 






3/2 


348.0593 


-0.2013 


-0 


1711 


-0 


0016 


-0 


0010 


347 


6843 


(10) 






5/2 


348.3697 


-0.2612 


-0 


1706 


-0 


0016 


-0 


0012 


347 


9351 


(10) 


1sCS)2s2pCP°) 


2po 


1/2 


348.0211 


-0.2212 


-0 


1506 


-0 


0005 


-0 


0001 


347 


6489 


(32) 







3/2 


348.1601 


-0.2430 


ls( o)2p ( U) 


2 n 


3/2 


349.5619 


—0.2113 






5/2 


349.7449 


-0.2810 




2p 


1/2 


349.8177 


-0.2749 






3/2 


350.4228 


-0.2681 


lsfS)2p2(i5) 




1/2 


351.7122 


-0.2338 


V (Z = 23), 1 - n/m = 


0.00001077, R = 


3.599 fm 






1/2 


-1152.2254 


0.3244 




2po 


1/2 


3.1007 


0.0257 






3/2 


3.8166 


-0.0034 






1/2 


376.9220 


-0.2932 


ls{^S)2s2p{^P'') 




1/2 


377.6474 


-0.2356 






3/2 


377.8797 


-0.2954 






5/2 


378.3762 


-0.3169 


1sCS)2s2pCP°) 


2po 


1/2 


380.2073 


-0.2783 






3/2 


380.6974 


-0.3053 


ls(25)2p2(3p) 


4p 


1/2 


381.2215 


-0.2314 






3/2 


381.5525 


-0.2312 






5/2 


381.9158 


-0.3008 


Isf5)2s2p(ip°) 


2po 


1/2 


381.4789 


-0.2569 






3/2 


381.6291 


-0.2776 


ls(25)2p2(i£)) 




3/2 


383.1346 


—0.2444 






5/2 


383.3790 


-0.3221 


lsfS)2p2(3p) 


2p 


1/2 


383.3748 


-0.3148 






3/2 


384.1218 


-0.3061 


\s{^S)2p'CS) 




1/2 


385.4416 


-0.2691 


Cr (Z = 24), 1 - n/m 


= 0.00001056, i? = 


= 3.642 fm 






1/2 


-1257.4911 


0.3705 


ls^2p 


2po 


1/2 


3.2560 


0.0296 






3/2 


4.1165 


-0.0038 


ls2s^ 




1/2 


411.6571 


-0.3348 


ls(^S)2s2p(3p°) 


4po 


1/2 


412.4160 


-0.2691 






3/2 


412.6880 


-0.3382 






5/2 


413.2899 


-0.3619 


ls(2S)2s2p(3p°) 


2po 


1/2 


415.1158 


-0.3143 






3/2 


415.6178 


-0.3494 


ls(2S)2p2(3p) 


4p 


1/2 


416.2168 


-0.2645 






3/2 


416.6254 


-0.2641 






5/2 


417.0454 


-0.3445 


Isf5)2s2p(ip'') 


2po 


1/2 


416.5140 


-0.2966 






3/2 


416.6968 


-0.3153 


ls('S)2p2(iD) 




3/2 


418.2864 


— 0.2809 






5/2 


418.6077 


-0.3669 


ls(2S)2p2(3p) 


2p 
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662 


8620 


(23) 


lsCS)2p\^P) 


2p 


1/2 


662.8170 


-0.6998 


-0 


4958 





0007 





0023 


661 


6244 


(22) 






3/2 


665.3388 


-0.6745 


-0 


4903 





0003 





0012 


664 


1755 


(23) 


lsCS)2p^CS) 




1/2 


666.7886 


-0.6239 


-0 


4856 


-0 


0010 


-0 


0018 


665 


6762 


(23) 


Ga (Z = 31), 1 - 


u/m 


= 0.00000796, R = 


3.997 fm 
























1/2 


-2127.5105 


0.8221 


1 


0703 





0001 


-0 


0003 


-2125 


6182 


(17) 


ls22p 


2po 


1/2 


4.3826 


0.0687 


-0 


0679 


-0 


0011 


-0 


0014 


4 


3809 


(23) 






3/2 


6.9654 


-0.0083 


-0 


0626 


-0 


0011 


-0 


0031 


6 


8903 


(25) 


Is2s2 




1/2 


698.9307 


-0.7403 


-0 


4288 


-0 


0002 





0002 


697 


7617 


(16) 


lsCS)2s2p(^P°) "P" 


1/2 


699.9339 


-0.5982 


-0 


4887 


-0 


0012 


-0 


0011 


698 


8448 


(14) 






3/2 


700.5811 


-0.7658 


-0 


4879 


-0 


0011 


-0 


0015 


699 


3250 


(20) 






5/2 


702.5455 


-0.8029 


-0 


4831 


-0 


0012 


-0 


0028 


701 


2554 


(11) 


U(^S)2s2p{^P°) 


1/2 


703.6159 


-0.6458 


-0 


4882 


-0 


0001 





0013 


702 


4831 


(218) 






3/2 


705.3997 


-0.7805 


-0 


4839 





0009 





0024 


704 


1386 


(13) 


ls(25)2p2(3p) 


4p 


1/2 


705.5012 


-0.5966 


-0 


5542 


-0 


0020 


-0 


0025 


704 


3458 


(22) 







3/2 


707.0412 


-0.5877 






5/2 


707.9046 


-0.7823 


1sCS)2s2pCP°) 


1/2 


706.5825 


-0.7024 






3/2 


706.9519 


-0.6755 


lsfS)2p2(i/)) 




3/2 


709.2370 


—0.6430 






5/2 


710.7834 


-0.7952 




2p 


1/2 


709.2091 


-0.7712 






3/2 


712.1194 


-0.7439 






1/2 


713.5797 


-0.6919 


Ge(Z = 32), l-n/m 


= 0.00000742, R = 


4.074 fm 


ls^2s 




1/2 


-2271.0881 


0.9075 






1/2 


4.5501 


0.0762 






3/2 


7.5083 


-0.0092 


ls2s^ 


2 c 


1/2 


740.3513 


-0.8168 


1sCS)2s2pCP°) 


4po 


1/2 


747.3914 


-0.6606 






3/2 


748.1030 


-0.8480 






5/2 


750.3834 


-0.8862 


ls(2S)2s2p(3p°) 


2po 


1/2 


751.2140 


-0.7065 






3/2 


753.3017 


-0.8612 


Isf5)2p2(3p) 


4p 


1/2 


753.2278 


-0.6610 






3/2 


755.0515 


-0.6492 






5/2 


755.9772 


-0.8662 


ls(2S')2s2p(^P°) 


2po 


1/2 


754.5352 


-0.7805 






3/2 


754.9348 


-0.7427 


ls(2S)2p2(i£)) 




3/2 


757.3212 


—0.7116 






5/2 


759.1781 


-0.8749 


lsfS)2p2(3p) 


2p 


1/2 


757.2496 


-0.8474 






3/2 


760.5881 


-0.8182 


15(^5)2^2(15) 




1/2 


762.0579 


-0.7639 


As (Z = 33), 1 - n/m = 


= 0.00000732, R = 


4.097 fm 


15^25 




1/2 


-2419.5616 


0.9986 


15^2^ 


2po 


1/2 


4.7194 


0.0843 






3/2 


8.0933 


-0.0101 


ls2s^ 




1/2 


795.3885 


-0.8983 


lsCS)2s2p(^P'') 


ipO 


1/2 


796.4664 


-0.7274 






3/2 


797.2441 


-0.9362 






5/2 


799.8803 


-0.9752 


lsCS)2s2p{^P°) 


2po 


1/2 


800.4294 


-0.7710 






3/2 


802.8605 


-0.9472 


lsCS)2p\^P) 


4p 


1/2 


802.5739 


-0.7302 






3/2 


804.7224 


-0.7149 






5/2 


805.7087 


-0.9559 


1sCS)2s2pCP°) 


2po 


1/2 


804.1468 


-0.8640 






3/2 


804.5775 


-0.8143 


ls(25)2p2(iD) 




3/2 


807.0554 


— 0.78OZ 






5/2 


809.2657 


-0.9609 


lsCS)2p^(^P) 


2p 


1/2 


806.9501 


-0.9282 






3/2 


810.7578 


-0.8977 


lsCS)2p^^S) 




1/2 


812.2371 


-0.8408 


Se(Z = 34), l-/L(/m = 


= 0.00000687, R = 


4.140 fm 


15^25 




1/2 


-2572.9594 


1.0958 


15^2^ 


2po 


1/2 


4.8906 


0.0929 






3/2 


8.7233 


-0.0112 


Is2s2 




1/2 


846.0488 


-0.9853 


lsCS)2s2p{^P°) ^P° 


1/2 


847.1672 


-0.7987 






3/2 


848.0123 


-1.0307 






5/2 


851.0472 


-1.0702 


lsCS)2s2p(^P'') 


2po 


1/2 


851.2706 


-0.8397 






3/2 


854.0877 


-1.0387 



-0 


5531 


-0 


0023 


-0 


0038 


705 


8943 


(24) 


-0 


5524 


-0 


0020 


-0 


0039 


706 


5640 


(24) 


-0 


4837 


-0 


0001 





0005 


705 


3967 


(83) 


-0 


4830 


-0 


0012 


-0 


0022 


705 


7900 


(131) 


-0 


5531 


-0 


0006 





0005 


708 


0407 


(30) 


-0 


5483 


-0 


0015 


-0 


0037 


709 


4348 


(26) 


-0 


5539 





0006 





0025 


707 


8871 


(24) 


-0 


5477 





0003 





0012 


710 


8294 


(26) 


-0 


5429 


-0 


0010 


-0 


0021 


712 


3418 


(24) 


1 


1906 





0001 


-0 


0003 


-2268 


9901 


(18) 


-0 


0758 


-0 


0011 


-0 


0015 


4 


5478 


(25) 


-0 


0698 


-0 


0011 


-0 


0035 


7 


4247 


(27) 


-0 


4768 


-0 


0002 





0003 


745 


0578 


(17) 


-0 


5439 


-0 


0012 


-0 


0012 


746 


1846 


(15) 


-0 


5430 


-0 


0010 


-0 


0016 


746 


7093 


(21) 


-0 


5376 


-0 


0012 


-0 


0033 


748 


9552 


(12) 


-0 


5435 


-0 


0001 





0014 


749 


9652 


(531) 


-0 


5384 





0008 





0027 


751 


9056 


(13) 


-0 


6169 


-0 


0020 


-0 


0027 


751 


9452 


(24) 


-0 


6156 


-0 


0022 


-0 


0043 


753 


7802 


(27) 


-0 


6151 


-0 


0019 


-0 


0044 


754 


4897 


(27) 


-0 


5382 


-0 


0001 





0007 


753 


2170 


(186) 


-0 


5377 


-0 


0012 


-0 


0025 


753 


6507 


(364) 


-0 


6157 


-0 


0006 





0006 


755 


9938 


(62) 


-0 


6101 


-0 


0015 


-0 


0043 


757 


6872 


(46) 


-0 


6165 





0006 





0027 


755 


7891 


(26) 


-0 


6095 





0003 





0013 


759 


1619 


(56) 


-0 


6046 


-0 


0010 


-0 


0025 


760 


6860 


(37) 


1 


3200 





0001 


-0 


0003 


-2417 


2431 


(19) 


-0 


0845 


-0 


0012 


-0 


0017 


4 


7164 


(27) 


-0 


0776 


-0 


0011 


-0 


0040 


8 


0004 


(29) 


-0 


5284 


-0 


0002 





0003 


793 


9619 


(19) 


-0 


6033 


-0 


0012 


-0 


0013 


795 


1332 


(20) 


-0 


6024 


-0 


0011 


-0 


0018 


795 


7026 


(29) 


-0 


5961 


-0 


0012 


-0 


0038 


798 


3039 


(13) 


-0 


6030 


-0 


0002 





0014 


799 


0566 


(550) 


-0 


5970 





0008 





0029 


801 


3201 


(14) 


-0 


6843 


-0 


0020 


-0 


0029 


801 


1544 


(25) 


-0 


6830 


-0 


0023 


-0 


0049 


803 


3173 


(37) 


-0 


6825 


-0 


0020 


-0 


0049 


804 


0634 


(67) 


-0 


5968 


-0 


0001 





0008 


802 


6867 


(218) 


-0 


5963 


-0 


0012 


-0 


0029 


803 


1627 


(488) 


-0 


6831 


-0 


0006 





0007 


805 


5863 


(998) 


-0 


6765 


-0 


0017 


-0 


0051 


807 


6216 


(1014 


-0 


6838 





0006 





0030 


805 


3416 


(28) 


-0 


6760 





0003 





0014 


809 


1858 


(195) 


-0 


6710 


-0 


0010 


-0 


0029 


810 


7213 


(106) 


1 


4588 





0001 


-0 


0002 


-2570 


4050 


(21) 


-0 


0937 


-0 


0012 


-0 


0018 


4 


8867 


(29) 


-0 


0860 


-0 


0011 


-0 


0046 


8 


6205 


(31) 


-0 


5837 


-0 


0002 





0003 


844 


4799 


(29) 


-0 


6669 


-0 


0012 


-0 


0014 


845 


6989 


(20) 


-0 


6661 


-0 


0010 


-0 


0019 


846 


3125 


(119) 


-0 


6589 


-0 


0012 


-0 


0043 


849 


3126 


(13) 


-0 


6667 


-0 


0002 





0015 


849 


7654 


(405) 


-0 


6597 





0008 





0032 


852 


3932 


(15) 





4p 


1/2 


853.5469 


-0.8045 


-0 


7565 


-0.0020 


-0.0031 


851.9808 


(27) 






3/2 


856.0652 


-0.7851 


-0 


7552 


-0.0023 


-0.0056 


854.5169 


(35) 






5/2 


857.1102 


-1.0518 


-0 


7547 


-0.0020 


-0.0056 


855.2962 


(234) 




) 


1/2 


855.4287 


-0.9531 


-0 


6596 


0.0000 


0.0010 


853.8170 


(198) 






3/2 


855.8914 


-0.8906 


-0 


6591 


-0.0012 


-0.0033 


854.3372 


(459) 






ji ^ 




u.ouo± 


n 

u 


I ooo 




n nnns 

u.uuuo 


S^fi Sfi'^^ 


, J 






5/2 


861.0919 


-1.0507 


-0 


7478 


-0.0017 


-0.0060 


859.2857 


(80) 


ls{^S)2p^('P) 


2p 


1/2 


858.3216 


-1.0140 


-0 


7560 


0.0005 


0.0032 


856.5553 


(30) 






3/2 


862.6461 


-0.9817 


-0 


7473 


0.0002 


0.0015 


860.9187 


(31) 


ls{^S)2p'^CS) 




1/2 


864.1319 


-0.9252 


-0 


7423 


-0.0010 


-0.0034 


862.4600 


(55) 


Br (Z = 35), 1 - 


/i/m : 


= 0.00000695, P = 


4.163 fm 














ls^2s 




1/2 - 


2731.3084 


1.1992 


1 


6072 


0.0001 


-0.0002 


-2728.5021 


(22) 


ls^2p 




1/2 


5.0639 


0.1021 


-0 


1038 


-0.0013 


-0.0020 


5.0589 


(31) 






3/2 


9.4015 


-0.0122 


-0 


0950 


-0.0012 


-0.0052 


9.2879 


(33) 


ls2s^ 




1/2 


898.3291 


-1.0778 


-0 


6428 


-0.0002 


0.0003 


896.6087 


(104) 


lsCS)2s2pi^P'' 


) '^p° 


1/2 


899.5015 


-0.8747 


-0 


7351 


-0.0013 


-0.0016 


897.8889 


(18) 






3/2 


900.4152 


-1.1315 


-0 


7342 


-0.0011 


-0.0021 


898.5463 


(73) 






5/2 


903.8953 


-1.1711 


-0 


7260 


-0.0013 


-0.0049 


901.9918 


(14) 


ls{^S)2s2p('^P° 


) 


1/2 


903.7457 


-0.9124 


-0 


7349 


-0.0003 


0.0015 


902.0997 


(321) 






3/2 


906.9946 


-1.1357 


-0 


7269 


0.0008 


0.0034 


905.1362 


(16) 


ls(2S)2p2(3p) 


4p 


1/2 


906.1546 


-0.8838 


-0 


8338 


-0.0021 


-0.0033 


904.4315 


(29) 






3/2 


909.0908 


-0.8599 


-0 


8325 


-0.0025 


-0.0063 


907.3896 


(33) 






5/2 


910.1947 


-1.1536 


-0 


8320 


-0.0021 


-0.0062 


908.2007 


(40) 


1sCS)2s2pCP° 


) 


1/2 


908.3921 


-1.0480 


-0 


7267 


0.0000 


0.0012 


906.6186 


(171) 






3/2 


908.8879 


-0.9717 


-0 


7263 


-0.0013 


-0.0038 


907.1849 


(396) 


1 '^"i9«2|'i n\ 


U 


3/2 


911.5806 


n QziRn 

— u.y^iuu 


n 




n nnnfi 


n nm n 


ono 8n9'^ 








5/2 


914.6402 


-1.1469 


-0 


8241 


-0.0019 


-0.0070 


912.6604 


(39) 


ls(2S')2p2(3p) 


2p 


1/2 


911.3776 


-1.1048 


-0 


8333 


0.0005 


0.0034 


909.4434 


(32) 






3/2 


916.2637 


-1.0713 


-0 


8236 


0.0002 


0.0016 


914.3705 


(34) 


\sQS)2p^('S) 




1/2 


917.7570 


-1.0135 


-0 


8185 


-0.0011 


-0.0039 


915.9200 


(53) 


Kr (Z = 36), 1 - 


li/m 


= 0.00000654, R = 


4.188 fm 


















1/2 - 


2894.6404 


1.3090 


1 


7657 


0.0001 


-0.0002 


-2891.5659 


(24) 


ls^2p 


2po 


1/2 


5.2393 


0.1120 


-0 


1143 


-0.0013 


-0.0022 


5.2334 


(33) 






3/2 


10.1312 


-0.0134 


-0 


1046 


-0.0012 


-0.0059 


10.0061 


(36) 


Is2s2 




1/2 


952.2849 


-1.1759 


-0 


7058 


-0.0002 


0.0003 


950.4035 


(25) 


\s(^S)2s2p(^P° 


) '^P" 


1/2 


953.4789 


-0.9554 


-0 


8078 


-0.0013 


-0.0017 


951.7127 


(19) 






3/2 


954.4617 


-1.2389 


-0 


8069 


-0.0011 


-0.0023 


952.4125 


(37) 






5/2 


958.4372 


-1.2784 


-0 


7977 


-0.0013 


-0.0056 


956.3541 


(15) 


\s(^S)2s2p(:^P° 


) 


1/2 


957.8644 


-0.9895 


-0 


8076 


-0.0003 


0.0016 


956.0686 


(268) 






3/2 


961.5943 


-1.2385 


-0 


7985 


0.0008 


0.0037 


959.5618 


(17) 


\s(^S)2p^(:'P) 


4p 


1/2 


960.4056 


-0.9684 


-0 


9164 


-0.0020 


-0.0036 


958.5153 


(30) 






3/2 


963.8120 


-0.9394 


-0 


9150 


-0.0025 


-0.0071 


961.9480 


(35) 






5/2 


964.9725 


-1.2619 


-0 


9146 


-0.0021 


-0.0070 


962.7869 


(36) 


\s(^S)2s2p{^P° 


) 


1/2 


963.0497 


-1.1487 


-0 


7983 


0.0000 


0.0015 


961.1041 


(147) 






3/2 


963.5798 


-1.0577 


-0 


7979 


-0.0013 


-0.0043 


961.7186 


(338) 


ls( o)2p ( U) 


2 n 
JJ 


3/2 


966.3745 


—1.0340 


— 


9151 


—0.0006 


0.0012 


964.4259 


36) 






5/2 


969.9362 


-1.2492 


-0 


9055 


-0.0019 


-0.0081 


967.7716 


(37) 


ls{^S)2p^CP) 


2p 


1/2 


966.1294 


-1.2011 


-0 


9158 


0.0004 


0.0037 


964.0167 


(34) 






3/2 


971.6277 


-1.1665 


-0 


9051 


0.0002 


0.0017 


969.5580 


(36) 


lsCS)2p\^S) 




1/2 


973.1282 


-1.1071 


-0 


8999 


-0.0011 


-0.0045 


971.1155 


(57) 



TABLE V: Theoretical wavelengths of the Is2l2l' —J- 13^21' transitions in lithlium-like ions, in A. 
Transitions are labelled as by Gabriel |39ll . 



Key 




Transition 






Z = 18 


Z = 19 


Z = 20 


Z = 21 


Z = 22 


Z = 23 


Z = 2i 


a 


ls(^S)2p2(3p) 




ls^2p 


2 po 

-^3/2 


3.985716(12) 


3.562881 (7) 


3.203642 (6) 


2.895849 (6) 


2.630130(5) 


2.399149 (5) 


2.197103 (4) 


b 


ls(^S)2p2(3p) 


'-P3/2 ^ 


ls^2p 


2 po 
-^1/2 


3.981655 (8) 


3.558756 (6) 


3.199456 (6) 


2.891608 (5) 


2.625838 (4) 


2.394811 (4) 


2.192723(4) 


c 


ls(^S)2p2(3p) 




ls^2p 


2po 
-^3/2 


3.989796 (22) 


3.567118 (7) 


3.208023 (7) 


2.900372 (6) 


2.634791 (6) 


2.403940 (5) 


2.202015 (5) 


d 


ls(^S)2p2(3p) 


'Pl/2 ^ 


ls^2p 


2po 


3.985727 (20) 


3.562983 (7) 


3.203826 (6) 


2.896118 (5) 


2.630484 (5) 


2.399585 (4) 


2.197616(4) 


e 


ls(2S)2p2(3p) 


-^5/2 


ls^2p 


2 po 
^3/2 


A.mnn (i2) 


3.586089 (7) 


3.223791 (6) 


2.913522(6) 


2.645790 (5) 


2.413162 (5) 


2.209764 (4) 



14 



/ 


ls(^S 


g 


IsCS 


h 


ls(^S 


i 


ls(^5 


3 


ls(^S 


k 


ls(^S 


I 


ls( S 


m 


ls( S 


n 




o 




P 


ls2s 


q 


1 ^^2 c- 


r 


IsCS 


s 


Is^S 


t 


IsCS 


u 


IsCS 


V 


isCs 



2p2(3p)''Pi/2 - 

2p=^(lD) 2^3/2 

2p2(i5)^Si/2 - 
2p2(i5)2Si/2 - 
■Si/2 ^ ls^2p2 

2s2p(^P°)' 
2s2p(3p°)2p, 

2s2p(ip<')^P°/2 

2s2p(3p'')*P3°/2 ■ 
2s2p(^P°) 



ls22p=^P°/2 
ls"2p"P°/2 
l,s22p2p°/2 
ls22p2po^^ 

> ls22p^P°^2 

> ls^2p^ 

> ls^2p^ 



iy2 

°3/2 



liS^2p^P; 



3/2 



ls^2p^P°/2 



1/2 ■ 



1/2 - 



ls^2s^5i/2 
ls^2s^Si/2 
ls^2s2Si/2 
ls^2s2Si/2 
ls^2s^Si/2 
U''2s^Si/2 



4.014740 
4.010620 
4.016579 
4.012454 
3.993809 
3.989864 
3.993942 
3.965728 
3.961707 
4.072729 
4.068489 
3.981328 
3.983465 
3.967689 
3.968543 
4.014997 
4.016105 



(11) 

(7) 

(12) 

(7) 

(61) 

(58) 

(59) 

(230) 

(229) 

(35) 

(33) 

(7) 

(15) 

(75) 

(68) 

(6) 

(6) 



3.588094 (7) 
3.583910(6) 
3.589979 (7) 
3.585791 (6) 
3.570066 (29) 
3.566231 (28) 
3.570374 (28) 
3.546235 (55) 
3.542148 (54) 
3.637985 (27) 
3.633684 (26) 
3.559148 (6) 
3.561332 (12) 
3.547584 (56) 
3.548452 (50) 
3.588101 (6) 
3.589219 (6) 



3.225781 (6) 
3.221537 (6) 
3.227718 (7) 
3.223469 (6) 
3.210076(16) 
3.206360(15) 
3.210564(16) 
3.189673 (27) 
3.185524(27) 
3.269158 (22) 
3.264800 (21) 
3.200463 (5) 
3.202698 (11) 
3.190554 (43) 
3.191430 (37) 
3.225603 (5) 
3.226725 (5) 



2.915487 (6) 
2.911188(5) 
2.917481 (6) 
2.913176(5) 
2.901656(10) 
2.898069 (9) 
2.902330 (10) 
2.884054(17) 
2.879847 (17) 
2.953565 (18) 
2.949154 (18) 
2.893146(5) 
2.895435 (9) 
2.884570 (34) 
2.885449 (29) 
2.915172(5) 
2.916294(5) 



2.647719 (5) 
2.643369 (5) 
2.649774 (6) 
2.645418 (5) 
2.635408 (7) 
2.631959(7) 
2.636271 (7) 
2.620117(13) 
2.615858 (13) 
2.681441 (15) 
2.676980 (15) 
2.627836(4) 
2.630186 (9) 
2.620349 (29) 
2.621228 (24) 
2.647301 (4) 
2.648416(4) 



2.415043 (5) 
2.410648 (4) 
2.417167 (5) 
2.412764 (4) 
2.403973 (6) 
2.400670 (6) 
2.405029 (6) 
2.390606 (10) 
2.386299 (10) 
2.445151 (13) 
2.440646 (12) 
2.396693 (257) 
2.399627 (9) 
2.390680 (42) 
2.391488 (20) 
2.414549 (4) 
2.415655 (4) 



2.211588(4) 
2.207150(4) 
2.213788 (5) 
2.209341 (4) 
2.201533 (5) 
2.198382 (5) 
2.202784 (5) 
2.189781 (8) 
2.185430 (8) 
2.238674(11) 
2.234127(11) 
2.195468(3) 
2.197951 (8) 
2.189613 (21) 
2.190472(17) 
2.211039(4) 
2.212130 (4) 



. 25 



26 



- 27 



28 



. 29 



: 30 



^ 31 



2p^CP)''P3/2 

2p2(-'*P)2p3/2 

2^='P)2Pl/2 
3p) 2p^ 



2p^ 
2p^ 



2p^CP)' 

2/('P)*P3/2 



^1/2 



^5/2 



3/2 



ls^2p 
ls^2p 
ls^2p 
ls^2p 
ls^2p 



2p2f P)*P, 
2p^(='P)*Pi 



-3/2 

^3/2 
P3/2 



2p^ 



/2 ■ 
°l/2 ■ 



l«^2p^ 
ls^2p2 



2p2(iD)^D5/2 
2p2(iD)2D3/2 
2p^CD)^Ds,2 
2p=^(iS)^Si/2 
2p2(iS)^Si/2 
■Si/2 ls^2p^P°^2 
'■Si, 2 ^ ls22p2p°/2 
2s2pC''P°)2p°^2 
2s2pC''P°)2p°^2 
2s2p(ip°)^P3°/2 
2s2p(ip°)2p°^2 ■ 
2s2p(''P°)^P3°/2 ■ 
2s2p(^ P°)*Pj 



ls^2p^P3°y2 

ls^2p^P°,2 
ls'-'2p^P°,2 

ls^2p2p3°/2 
ls=^2p=^P°/2 



ls^2s^Si/2 
ls^2s^Si/2 
ls^2s\Si/2 
ls^2s^S 



l/2_ 



1/2 

ls''2s^5i/2 
ls^2a^Si/2 



2.019352(4) 
2.014933 (3) 
2.024374 (4) 
2.019934(4) 
2.030897 (4) 
2.032655 (4) 
2.028178 (4) 
2.034937 (4) 
2.030450 (4) 
2.023433 (4) 
2.020438 (4) 
2.024881 (5) 
2.013050(7) 
2.008658 (6) 
2.057194 (9) 
2.052608 (9) 
2.017983 (3) 
2.020539 (7) 
2.012747 (19) 
2.013592 (15) 
2.032068 (3) 
2.033141 (3) 



1.862150(4) 
1.857696 (3) 
1.867272(4) 
1.862793 (3) 
1.872775 (4) 
1.874458 (4) 
1.869945 (3) 
1.876829(4) 
1.872304(4) 
1.865920(4) 
1.863083 (4) 
1.867564(4) 
1.856704(6) 
1.852275 (5) 
1.896835 (8) 
1.892213 (8) 
1.861013 (3) 
1.863646 (7) 
1.856299 (17) 
1.857127 (13) 
1.873847 (3) 
1.874897 (3) 



1.722447 (3) 
1.717960 (3) 
1.727654(4) 
1.723139(3) 
1.732309(4) 
1.733912(4) 
1.729364(3) 
1.736377 (4) 
1.731816(3) 
1.725934(4) 
1.723258 (4) 
1.727773 (4) 
1.717723 (5) 
1.713259 (5) 
1.754442 (7) 
1.749786 (7) 
1.721512(3) 
1 .724224 (6) 
1.717241 (15) 
1.718049(11) 
1.733287 (3) 
1.734312 (3) 



1.597740 (3) 
1.593221 (2) 
1.603022(4) 
1.598473(3) 
1.606968 (3) 
1.608485 (3) 
1.603905 (3) 
1.611051 (4) 
1.606457 (3) 
1.600968 (3) 
1.598451 (3) 
1.603000(4) 
1.593626(4) 
1.589130(4) 
1.627428 (7) 
1.622741 (7) 
1.596979 (3) 
1.599773 (5) 
1.593087 (12) 
1.593874(9) 
1.607856(3) 
1.608852(3) 



1.485954 (3) 
1.481407 (2) 
1.491298 (3) 
1.486718 (3) 
1.494655 (3) 
1.496082 (3) 
1.491473 (3) 
1.498754 (4) 
1.494129(3) 
1.488945 (3) 
1.486586 (3) 
1.491165 (4) 
1.482358 (4) 
1.477833 (3) 
1.513656(6) 
1.508938 (6) 
1.485345 (2) 
1.488224(3) 
1.481783 (8) 
1.482547 (5) 
1.495458 (3) 
1.496423 (3) 



1.385365 (3) 
1.380792 (2) 
1.390759 (3) 
1.386150 (3) 
1.393626 (3) 
1.394963 (3) 
1.390326 (2) 
1.397745 (3) 
1.393089 (3) 
1.388137 (3) 
1.385932 (3) 
1 .390540 (4) 
1.382212 (3) 
1.377659 (3) 
1.411347 (6) 
1.406601 (5) 
1.384886 (2) 
1.387849 (9) 
1.381603 (3) 
1.382350 (3) 
1.394349 (3) 
1.395280 (3) 



1.294526(3) 
1.289927(2) 
1.299959 (3) 
1.295322 (2) 
1.302417 (3) 
1.303665 (3) 
1.299002 (2) 
1.306559 (3) 
1.301875 (3) 
1.297095 (3) 
1.295039(4) 
1 .299674 (4) 
1.291750(3) 
1.287171 (2) 
1.319011 (5) 
1.314238 (5) 
1.294159 (2) 
1 .297208 (40) 
1.291131 (24) 
1.291850(15) 
1.303067 (4) 
1.303962 (3) 



Transition 



Z = 32 



Z = 33 



Z = 34 



Z = 35 



Z = 36 



2p- 
2p^ 
2p^ 
2p^ 
2p^ 
2p^ 
2p^ 
2p^ 
2p2 
2p= 
2p= 
2p2 
2p2 

y 

'Sl/2 



)'P3/2 ^ l.S^2p"P° 

,2^1s^2p^P° 
,2 ls^2p^P° 



f^5/2 ■ 
P3/2 



r^3/2 ■ 



ls^2p^ 
ls^2p^ 
ls^2p^ 



3/2 

°l°/2 



1/2 - 

'0)=D3/2 

^D)^D3,2 

^S) ^Si/2 - 
ls^2p^ 



ls^2p2p°/, 

i. ls^2p2p3''/2 

> ls=2p2p°/2 



ls''2p^Pj 



ls^2p^ 
lis^2p^ 



3/2 
P3/2 



1/2 



Si 



/2 



Is 



2p^ 
PI 

2s2pC'P°)^P°/2 ■ 
2s2p(ip°)2p3°/2 
2s2p(ip°)2p°^2 
2s2p(3p°)*P°^, 

2a2p(''P°)^P°^2 



ls''2s^Si/2 
ls^2s^Si/2 
ls^2s^5i/2 
ls22s^5i/2 
ls^2s^Si/2 
ls^2s^Si/2 



1.212215 (8) 
1.207593 (8) 
1.217678 (3) 
1.213015(2) 
1.219796 (3) 
1.220956 (3) 
1.216268(2) 
1.223965 (3) 
1.219254 (2) 
1.214598 (7) 
1.212685 (9) 
1.217345 (10) 
1.209762 (5) 
1.205159 (5) 
1.235393 (5) 
1.230594(5) 
1.211943 (2) 
1.215079 (86) 
1.209137 (58) 
1.209833 (30) 
1.220377 (3) 
1.221235 (2) 



1.137399 (27) 
1.132755 (27) 
1.142882 (3) 
1.138194(2) 
1.144717(9) 
1.145791 (4) 
1.141079(4) 
1.148916(3) 
1.144178 (2) 
1.139624(144) 
1.137847 (142) 
1.142532 (143) 
1.135223 (15) 
1.130597 (15) 
1.159430(5) 
1.154605(4) 
1.137207 (2) 
1.140429 (79) 
1.134598 (69) 
1.135271 (31) 
1.145236(4) 
1.146056(3) 



1.069188 (2) 
1.064525 (2) 
1.074690(3) 
1.069978 (2) 
1.076288 (30) 
1.077280 (3) 
1.072546 (3) 
1.080519(3) 
1.075757 (2) 
1.071240(9) 
1.069591 (6) 
1.074300(6) 
1.067258 (6) 
1.062611 (6) 
1.090216(5) 
1.085367 (5) 
1.069069 (2) 
1.072375 (51) 
1.066636 (57) 
1.067286 (25) 
1.076750(15) 
1.077531 (3) 



1.006833 (2) 
1.002151 (1) 
1.012344(2) 
1.007610(2) 
1.013744(4) 
1.014659 (3) 
1.009904(2) 
1.018012(3) 
1.013225 (2) 
1.008739(3) 
1.007211 (2) 
1.011941 (3) 
1.005112(5) 
1.000446(5) 
1.026987 (13) 
1.022116(12) 
1.006773 (2) 
1.010162 (36) 
1.004500(44) 
1.005127 (19) 
1.014157 (8) 
1.014900(2) 



0.949680 (2) 
0.944980(1) 
0.955196(2) 
0.950441 (2) 
0.956429 (2) 
0.957272 (2) 
0.952496 (2) 
0.960736 (3) 
0.955926 (2) 
0.951451 (2) 
0.950036 (2) 
0.954786 (2) 
0.948141 (5) 
0.943456 (5) 
0.969023 (4) 
0.964130(4) 
0.949670 (2) 
0.953140(27) 
0.947540 (33) 
0.948146(15) 
0.956799 (4) 
0.957502 (2) 
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TABLE VI: Energy levels of Ar XVI {Z = 18), relative to the ground state. Comparison is made with the recent compilation of theoretical 
and experimental world data by Saloman |40]. Units are Rydbergs. 



State 


J 


Present work 


Saloman [40] 




ill 


2.3417(6) 


2.342 193 (36) 




3/2 


2.5749 (8) 


2.575 268 (36) 




1/2 
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226.9032 (4) 
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3/2 
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1/2 
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229.467(12) 
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TABLE VII: Energy levels of Fe XXIV {Z = 26), relative to the ground state. Comparison is made with the compilation of theoretical and 
experimental world data by Shirai et al. | 41] and with recent calculation by Safronova and Safronova |21J. Units are Rydbergs. 



State 


J 


Present work 


Shirai [4r] 


Safronova [211 


Is22p2p° 


ill 
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3/2 


4.7442 (16) 
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4.749 
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490.919 




5/2 


491.3308(16) 
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TABLE VIII: Energy levels of Rr XXXIV {Z — 36), relative to the ground state. Comparison is made with the compilation of theoretical and 
experimental world data by Saloman |42] and with recent calculation by Safronova and Safronova |2I]. Units are Rydbergs. 
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